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We present a calculation of the SU(3)-breaking corrections to the hyperon transition vector form 
factors to Cl(p^) in heavy baryon chiral perturbation theory with finite-range regularisation. Both 
octet and decuplet degrees of freedom are included. We formulate a chiral expansion at the kinematic 
point = —{Mbi — which can be conveniently accessed in lattice QCD. The two unknown 

low-energy constants at this point are constrained by lattice QCD simulation results for the —>■ n 
and —>■ transition form factors. Hence we determine lattice-informed values of fi at the 

physical point. This work constitutes progress towards the precise determination of |Kis| from 
hyperon semileptonic decays. 

PACS numbers: 12.38.Gc, 14.20.Jn, 12.39.Fe 

Keywords: Hyperon vector form factors, Lattice QCD, Chiral Extrapolation, Semi-leptonic decays, CKM 


I. INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa (CKM) quark¬ 
mixing matrix elements are fundamental parameters of 
the Standard Model. The precise determination of these 
quantities, which encode the flavor structure of the quark 
sector, is thus of great importance [1]. In particular, a 
stringent test of CKM unitarity [5] is given by the first- 
row relation |t4dP + |I4sP + |I4bP = 1- The matrix 
element | Vus \ has traditionally been extracted from mea¬ 
surements of kaon semileptonic and leptonic decays and 
the hadronic decays of tau leptons. These determinations 
have (until recently [MS]) been in slight tension |^. A 
further independent extraction of |I4s| can be performed 
based on hyperon beta decay. In particular, the prod¬ 
uct \Vusfi{Q^ = 0)1 can be extracted from experiment at 
the percent-level |7]; precise calculations of the hadronic 
corrections to the vector form factors fi{Q'^ = 0) are 
therefore of great interest. 

While the Ademollo-Gatto theorem | 8 ] protects the 
vector form factors from leading SU(3)-symmetry- 
breaking corrections generated by the mass difference 
of the strange and nonstrange quarks, knowledge of the 
second-order breaking corrections to fi{Q^ = 0 ) is cru¬ 
cial to obtain a precise value of |I4s| [3121 • Estimates of 
these corrections have been performed based on quark 
models [TOl E], l/Nc expansions [12], chiral effective 
field theory UMI] and quenched and unquenched lattice 
QCD [THHID]- An enduring puzzle has been that the sign 
of the SU(3) breaking corrections determined in lattice 
QCD (at this stage away from the physical pseudoscalar 
masses) and quark models is, in general, opposite to that 
determined from relativistic and heavy baryon chiral per¬ 
turbation theory and 1 /iVc expansions. 

Using finite-range regularised chiral perturbation the¬ 


ory, which appears to offer markedly improved con¬ 
vergence properties of the (traditionally poorly conver¬ 
gent [211 [22] ) SU (3) chiral expansion, we revisit the ef¬ 
fective held theory estimates of the SU(3) breaking in 
fi- To enable us to use 2-1-1-flavor lattice QCD to inform 
this determination, we choose to investigate the form fac¬ 
tors at = —{Mbi — Mb^Y (i- 6-5 9 = 0 in our lattice 
simulations with hxed zero sink momentum) instead of 
at = 0 as is standard. This choice avoids the need 
to model the Q^-dependence of the lattice simulation re¬ 
sults for fi and correct to = 0 before chiral and Hnite- 
volume corrections have been applied. This small shift 
in would be highly model-dependent with our simu¬ 
lation results; for example, after using linear and dipole 
interpolations to = 0 the results differ by more than 
the statistical uncertainty. 

After extrapolation to the physical regime we use simi¬ 
lar parameterisations to those used to ht experimental re¬ 
sults to extrapolate to = 0 , allowing comparison with 
other works. The sign of the SU(3)-symmetry-breaking 
corrections that we find is consistent with the results of 
quark models. Furthermore, while we have lattice simu¬ 
lation results for only the —)■ n and — )■ E+ tran¬ 

sitions, our chiral extrapolations also allow us to extract 
(lattice-informed) results for the SU(3)-breaking in the 
A ^ p and —>■ A transitions. 


II. CHIRAL PERTURBATION THEORY 

We use a formalism based on heavy-baryon chiral per¬ 
turbation theory, with hnite-range regularisation, to per¬ 
form chiral and infinite-volume extrapolations of lattice 
simulation results for fi. Our chiral expansion is different 
from those used previously HMZj in that we consider the 


2 


form factor at Q^ax = (corresponding to 

(f = 0 in our lattice simulations), rather than = 0. In 
addition to removing the model dependence of a —>■ 0 
interpolation, a feature of this approach is that several 
free low-energy constants appear in the chiral expansion 
(while fi{Q^ = 0) is protected from such coefficients to 
this order by Ademollo-Gatto 0). These free constants 
are fit to the lattice QCD simulation results. In practice, 
fitting these terms to the lattice simulation results may 
also absorb the effect of potentially significant higher- 
order terms neglected in the conventional chiral expan¬ 
sion. 

The structure of the master formula for our chiral ex¬ 
trapolation to 0{p^) is 

= [Mb-Mb'Y) 

=/fu(3) + _|_ (Tadpole Loops) 

-I- (Octet Loops) -I- (Decuplet Loops) 
-I- (I/Mq corrections). (1) 

The terms SU(3)-symmetric values of the 

vector form factors. At this order, the only undetermined 
chiral coefficients (which we will fit to the lattice simu¬ 
lation results) appear in the term proportional to Q^. 
These are the same coefficients which appear in expres¬ 
sions for the octet baryon electric radii (labelled bo and 
6^ in Ref. [23, for example). We discuss each remaining 
term of Eq. Q in turn. 


A. Terms linear in 


not contribute at = q and for this reason have not 
been considered by other authors. The following sections 
give expressions for the contribution of each diagram to 

Eq. Q. 


1. Octet, Decuplet and Tadpole loops 


For a transition Bi —>■ i? 2 , the wavefunction renormal¬ 
isation diagrams (Figs. 1(a) and|l(b)[) give contributions 


+Cb2B’4>^wr{B2, B')) , 

(4) 

{ I for octet baryon intermediate states 
1 for decuplet baryon intermediate states, 

(5) 

and a sum is implied over the repeated indices B' (repre¬ 
senting an octet or decuplet baryon intermediate state) 
and (j) (representing the intermediate meson). Here 


.^WR = 

where 
Xr = 


The term in Eq. Q which is linear in (i.e., for the 
application to our lattice simulations with g = 0 it is pro¬ 
portional to {Mb^ — Mb^Y): comes from the Lagrangian 
piece 

C = -ev>^{D''F+) [2bDTrB{X,B} + 2bFTTB[X,B]] , 

( 2 ) 

which is equivalent to Eq. (4) in Ref. [53] (where we have 
dropped the term proportional to TrA which vanishes 
here). Here A is the strangeness-changing matrix 

0 0 l\ 

000. (3) 

0 0 oy 

The two undetermined coefficients bo and bp, which ap¬ 
pear in the coefficients CbiB 2 which are tabulated in ap¬ 
pendix]^ will be fit to the lattice simulation results. We 
note that in Ref. |23| . the identically-named parameters 
are interpreted as chiral limit form factors at fixed values 
of (that is, they may in principle vary with Q^). In 
the present work, these terms may be interpreted as gen¬ 
uine chiral low-energy constants, i.e., the = 0 limit of 
those in Ref. [53] . 

B. Loop diagram contributions 


The meson loop diagrams included in this calculation 
are shown in Fig. a Note that Figs. |1(1)| and |l(m)| do 


UJ = 




(7) 


The baryon-baryon-meson coupling constants Cbb'^i 
are standard expressions (see, for example. Ref. |53| 1 
given in terms of the constants D, F and C. The tree- 
level vector current transition coefficients ff^BiB 2 ) 
the same as those in Eq. Q and are tabulated in ap¬ 
pendix For our numerical results the masses of the 
relevant baryons B' (at any values of TOtt and ttik ), which 
appear in the term involving octet baryon intermediate 
states (Fig. 1(a)), are taken from the fits to lattice simula¬ 
tion results for the octet baryon masses given in Ref. . 
For loops involving decuplet baryon intermediate states 
(Fig. |l(b)] ), the decuplet baryon masses (i.e., Mb> in the 
above equation) are set to the value 


Ma = 


■A 


phys. 


( 8 ) 


CTT/O'^ 

where ^ is the SU(3)-symmetric mass of the baryon 
octet at the particular meson masses of relevance, and 
Aphys. = 0.292 GeV is fixed. As discussed in detail 
later, this choice for A does not affect our final results; 
in fact, we find entirely consistent results for /i for each 
transition when contributions from intermediate decuplet 
states are omitted entirely. 
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FIG. 1. Loop diagrams included in the present calculation. Solid squares denote the usual strong-interaction meson-baryon 
vertices. The open circles and crosses represent 1 /Mq corrections (see Eq. ®) and insertions of the vector current, respectively. 
Single, double, and dashed lines denote octet baryons, decuplet baryons and mesons, respectively. 


Figures |l(c)| and |l(d)[ with vector current insertions into the intermediate baryons in the loop diagrams, generate 
contributions which can be expressed as 


Lr = 


B\ B"), 


(9) 


where 


Ib{B^,B2,B\B") 


1 r^kf _ 

3 Jo ~ {^B2 ~ Mb")){uj — {Mbi 


Mb')) 


Xb is as above, sums are again implied over the indices B', B" and (p, and the baryon masses are set as before in our 
numerical work. 


The contributions from the loops shown in 


Figs. 1(e) and |l(f)| 


can be written as 


Lm = 


Xr S2 , s'), 




( 10 ) 
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where 




4 + (Mb, - Mb,)) _ 

3 Jo yW 2 (wi +W 2 + (-^Bi — Mb 2)){^2 — Wi + {Mb^ — Mb 2)){^2 — {Mb 2 — Mb')) 

_ k*{2uji - {Mb2 - Mb2)) _\ 

Wi(wi +LU2 — (Mbi — MB2))i^l — UJ2 — {Mbi — Mb 2)){^^1 — {^Bi — Mb')) j 

( 11 ) 


and wi and W 2 denote the expression given in Eq. 0 
with the masses of the intermediate mesons and 
replacing m, respectively. Again, a sum is implied over 

CTT/'O'V 

s', (j)i and (j) 2 - The fi(^^\J, 2 ) tabulated in appendix 


The tadpole loop diagram shown in Fig. |l(g)| does not 
have any baryon intermediate states, so our expression 
reduces to that of Villadoro m if dimensional regulari- 
sation is used. We can write the contribution from this 
loop as 


Ex 


1 1 

8^2 


CbiB2<i>v It, 



where 


It 



(13) 


The coefficients CbxB 24 >i 4 > 2 V are given in ap pend ix 
Finally, the tadpole loop shown in Fig. |l(h)| gives a 
contribution 


1 1 

8^2 


^ .SU(3) 

'-'BiB24>i4>2 /i(0i02) 


.^TM, 


(14) 


where 


4tm = 2 



{2uji + (Mbi — Mbj))^ 

UJi(u!i + UJ2 + (Mbi — Mb2)){^1 — W2 + (ATbj — MB2)) 


(2uJ2 — (Mbi — Mbj))^ 

UJ2 (wi + UJ2 — {Mbi — Mb2 )) (wi — a;2 + {Mb^ 



(15) 


and all notation is as defined previously. The coefficients CbiB 2 (Pi 4>2 are again standard coefficients which may be 
found, for example, in Ref. [24] . 


2. 1 /Mo corrections 


The terms labelled as Oil/Me)) corrections are rela¬ 
tivistic corrections to the heavy-baryon formalism. They 
are generated by the Lagrangian pieces 


=^(B [{v^^dn)^ - d^d,] B) (16a) 

_ _ 

+ —D{d^BS^{A,2, B} - BS>^{A,, 

(16b) 

'1'})^ _ _ 

+ —F{d))BS^[A,2,B] - BS>^[A,,d)2B]). 

Mq 

(16c) 


The contribution from Eq. (16a I corresponds to (I/Mq) 
corrections to the baryon propagator, which we expand 


as 

i i 

k ■ V + ie k ■ V + ie 

i I .k'^ — {k ■ w)^\ i 
k ■ V + ie \ 2Mo J k ■ v + ie' 

(17) 

The remaining terms correspond to corrections to the 
strong-interaction vertices. Here Mq is the heavy-baryon 
mass scale, taken to be the octet baryon mass in the chi¬ 
ral limit. We use the value of Mq obtained in Ref. [53] 
from a chiral extrapolation of octet baryon masses. The 
difference between choosing the chiral-limit value or any 
physical baryon mass for Mq is a higher-order effect. In 
practice we do find that the shift in our results from 
making such a change of convention is small. As de¬ 
tailed later, we add a systematic accounting for this to 
the quoted uncertainties in our numerical work. 

Our expressions for the 0{1/Mq) correction terms 
agree with Villadoro m in the limit of no mass difference 
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between the baryons, no momentum transfer, and using 
dimensional regularisation to evaluate the loop integral 
expressions. We note again that the additional diagrams 
in Figs. 1(1) and l(m) contribute to our expressions but 
vanish at = 0. We do not consider decuplet inter¬ 
mediate states here. We write out each contribution in 
turn. 


The wavefunction renormalisation contribution shown 
in Fig. l(i) takes exactly the same form as Eq. Q, under 


the replacement Jwr —t Iwrm, where 
4 


IwRuiB, B') = — 


3Mt 


dk 




0 Jo \uj{uj — {Mb — Mbi)Y 

( 18 ) 

The total contribution from the other wavefunction 
renormalisation terms, depicted in Figs. l(j) and |l(k)[ 
also takes exactly the same form as Eq. Q, under the 
replacement Iwr —t .Iwrmv, where 


Iwkmy{B, B') = — 


3M, 


dk 


0 Jo 


{uj — {Mb — Mb')Y 
(19) 


The total contribution from the diagrams shown in Figs. mi and |l(m)| can be expressed in the form of Eq. ( |10[ ) , 
under the replacement /m —t /mmV; where 


Iumw{Bi,B2tB') = 

4 1 


3 4Mo Jo 


dk\ 


k\2uJ2 + {MB,-MB^)f 


UJ2{ui2 + UJi + {Mbi — Mb2)){(^2 — CJi -f {Mbi — Mb2)){(^2 — {Mb2 — Mb')) 

_ k*{2uJi - {Mb, - Mb2))^ _\ 

wi ( a;i + a ;2 — {Mb, — Mb2)){^i — ui2 — {Mb, — Mb2)){^i — {Mb, — Mb')) I 


( 20 ) 


Similarly, the contribution from Fig. l(n) can be expressed in the form of Eq. (101, under the replacement Im —t /mm, 
where 


^mm{Bi, B 2 , B') = 

4 1 _ k^2LU2 + {Mb, - MB 2 )) _ 

3 2Mo Jo y W2(a;2 -I- wi -|- {Mb, — Mb2)){>^2 — uji + {Mb, — Mb2)){^'J2 — {Mb2 — Mb'))"^ 

^ _ k^{2uji - {Mb, - MB 2 )) _ 

wi(a;i + UJ2 — {Mb, — Mb2)){^i —002 — {Mb, — MB2)){i^i — {Mb, — Mb'))^ 

The total contribution from Figs. mu and |l(p)] can be expressed in the same form as Eq. ([^, under the replacement 
Ib Ibm, where 



2 f 


k’ 


Lu{uj — {Mb, — Mb'))‘^{uj — {Mb2 — Mb")) 
k^ 

u}{oj — {Mb, — Mb')){uj — {Mb2 — Mb"))"^ 


( 22 ) 


The contribution from Figs. |l(q)| and |l(r) can also be expressed in the same form as Eq. ([^, under the replacement 
Ib -t Ibmy, where 




k^ 


{lo — {Mb, — Mb')){uj — {Mb2 — Mb")) 
Finally, the contribution from Figs. |l(s)1and|l(t)]can be expressed as 


T OCt _ 

— 


V - o_2 2 


8^V^ 


2 i?') - C'B^B'0C'3,g^^y/M(Si,O,i?',O) 


(23) 


(24) 


3. Master formula expressions for the loop diagram contributions: 

The master formula for the chiral extrapolation, 

Eq. 0, can now be re-written in terms of the explicit 
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= {Mb - MB>f) 

+ Cb,b,Q^ + Lt 

I T oct I T oct I T dec I r dec 
-^WR -^WR -^M 

I I r oct I r oct I r oct i r oct i roct 

-h -hi^wRM -^WRMV -^BM -^EMV -^V * 

(25) 


in the 1 /Mq correction diagrams (which account for rela¬ 
tivistic corrections to the heavy-baryon approach), which 
was not done previously. Our expressions match those 
from Villadoro m, which is the most closely related cal¬ 
culation, in the relevant limits. 

III. LATTICE SIMULATION 


In practice, we evaluate each loop integral expression 
within the finite-range regularisation (FRR) scheme. 
This choice of regularization procedure is discussed in 
detail in Refs. [2HVl28| . In short, the inclusion of a fi¬ 
nite cutoff into the loop integrands effectively resums the 
chiral expansion in a way that suppresses the loop contri¬ 
butions at large meson masses. This enforces the phys¬ 
ical expectation, based on the finite size of the baryon, 
that meson emission and absorption processes are sup¬ 
pressed for large momenta [53]. For the case of the octet 
baryon masses, FRR appears to offer markedly improved 
convergence properties of the (traditionally poorly con¬ 
vergent [ 2 TJI 22 I) SU(3) chiral expansion, and this scheme 
consistently provides robust fits to lattice data at lead¬ 
ing or next-to-leading order. Nevertheless, one could 
calculate the size of higher order corrections to confirm 
that these contributions are small as expected. Changing 
to dimensionally regularized integral expressions requires 
simple substitutions; details are given in [55] . 

For our numerical analysis we choose a dipole reg¬ 
ulator u{k) = ^ ^ with a regulator mass A = 

0.8 ± 0.2 GeV. This form is suggested by a comparison 
of the nucleon’s axial and induced pseudoscalar form fac¬ 
tors and the choice of A is supported by a lattice analysis 
of nucleon magnetic moments EHIISI]. We note that dif¬ 
ferent regulator forms, for example monopole, Gaussian 
or sharp cutoff yield extrapolated results which are con¬ 
sistent within the quoted uncertainties. 

The most recent comparable effective held theory stud¬ 
ies are those by Villadoro [T3], which is a heavy-baryon 
calculation (with dimensional regularisation) with both 
decuplet intermediate states and relativistic (I/Mq) cor¬ 
rections considered, Lacour which is a covariant 

calculation that neglects the decuplet, and Geng m, 
which is a covariant calculation with decuplet degrees of 
freedom. The primary difference between these works 
and ours, other than our non-zero Q^, is that we have 
more carefully accounted for the meson-mass dependence 
of the mass-splittings among the octet baryons. While 
other authors have either hxed the mass differences be¬ 
tween members of the baryon octet (which appear in 
the loop integral expressions) absolutely, or varied them 
linearly, we have used the chiral extrapolation of octet 
baryon masses which was presented in Ref. [5^ to give 
these mass differences as a function of meson masses. 
Furthermore, while Villadoro Taylor-expands in the mass 
differences between the octet baryons, we include them 
explicitly (i.e., resumming Villadoro’s expansion to all or¬ 
ders). We also include the octet baryon mass splittings 


We apply the formalism developed in the previous sec¬ 
tions to new lattice simulation results for the E —V 
and S —>■ E transition vector form factors. Preliminary 
results appeared in Refs. [351 [33] . 

We use gauge field configurations with 2-1-1 flavours of 
nonperturbatively 0(a)-improved Wilson fermions. The 
clover action consists of the tree-level Symanzik improved 
gluon action together with a mild ‘stont’ smeared fermion 
action [33]. We use a single lattice volume, L^T = 32^ x 
64, with f3 = 5.5, corresponding to the lattice scale a = 
0.074(2) fm (set using various singlet quantities [33U3H] V 
Details are given in Table |lj The simulations correspond 
to two distinct fixed-singlet-mass trajectories in rriT^-mx 
space. 

The vector form factor fi is defined in terms of the ma¬ 
trix elements of the charged strangeness-changing (s —>■ 
u) weak vector current = u'y^s. In Euclidean space, 
defining q'^ = —Q"^, 


{B2\VM = B2{P2) 


7m/i ( 9 ^ 




Mbi -I- Mb. 


■/2 




Mbi + Mb. 


■h ( 9 ^ 


Bi{pi). 

(26) 


While here we are interested in the vector form factor fi 
only, what is accessible on the lattice at = (-^Si ~ 
Msj)^ is in fact the ‘scalar form factor’: 

/o(.-) = /i(.-)+ ^. /3(.-). (27) 

Bl ±52 


This quantity can be obtained with high precision from 
the ratio 


R{t,t') 


Gf^^(t',t;0,0)Gf^^(F,t;0,0) 

Gf^i(t',t;0,0)Gf^^(t',t;0,0) 


(28) 


which tends to |/o(9max)|^ in th® limit t,{t' — t) —)■ 00 . 
Here, for example, Gf^^^{t',t]0,0) is the zero three- 
momentum lattice three-point function of the fourth 
component of the vector current inserted at time 
t between the source baryon Bi located at time t = 0 
and the sink baryon B 2 at time t'. In the SU(3) flavour 
symmetric limit, R{t',t) = I; any deviations from unity 
are purely due to symmetry-breaking effects. 

Away from = {Mb^ — Mb.^)'^^ it is possible to de¬ 
termine both fi and /a independently [151150] . The ratio 









7 


Ko Ks m.n- (MeV) rriK (MeV) rrinL 

/f" 

q^ (GeV^) lattice value 

q^ (GeV^) lattice value 

1 0.120900 0.121040 0.120620 360 505 4.3 

2 0.121095 0.120512 310 520 3.7 

0.0017 -1.000(3) 

0.0042 -0.999(7) 

0.0006 0.990(3) 

0.0005 0.985(5) 

3 0.120950 0.121040 0.120770 330 435 4.0 

0.0006 -1.002(4) 

0.0002 0.994(3) 


TABLE I. Details of the lattice simulation parameters and raw lattice simulation results for the hyperon vector transition form 
factors at g = 0 with fixed zero sink momentum. The parameter kq denotes the value of ki = Ks at the SU(3)-symmetric point, 
/? = 5.5 corresponding to a = 0.074(2) fm, and x T = 32^ x 64. 



-q^ (GeV^) 



FIG. 2. Ratio of /s to fi for the two transitions for which 
we produce lattice simulation results. The three sets of data 
(brown, red, purple) on each plot correspond to the three 
sets of pseudoscalar masses (labelled 1, 2, 3) in Table The 
vertical green band shows the range of values for which we 
wish to estimate / 3 // 1 , while the horizontal blue band shows 
our bounds on this ratio at the values of interest based on 
the linear extrapolations of the simulation results which are 
shown as brown, red and purple bands. 


of these two quantities, at the pseudoscalar masses of our 
simulations, is shown in Fig.[^ At the relevant values of 
9max = (-^Bi ~ indicated by the green vertical 

band on the figure, it is clear that I/3//1I < {0.15,0.2} 
(where the two numbers are for the S —E and Y, ^ N 
transitions respectively) is a conservative bound. This is 
supported by both quenched lattice simulations |T8l [19] 
and quark models m which find small values for that 
ratio over a range of values of q^. Taking these bounds, 
scaled by the appropriate kinematic factors, the contri¬ 


bution from /a to /o is then for each set of simulation 
pseudoscalar masses negligible compared with the sta¬ 
tistical uncertainties of the calculation. That is, to the 
statistical precision of our simulation, /q and fi are the 
same, so we take our simulation value of /o as the result 
for the vector form factor. 

We access the strangeness-changing transitions involv¬ 
ing only the outer-ring octet baryons, i.e., E —> fV and 
S —)■ E. Simulation results are given in Table We 
note that the simulations suffer no systematic uncertain¬ 
ties from omitted disconnected loops because these terms 
cannot contribute to the transitions. 


IV. FITS TO THE LATTICE RESULTS 

Before performing a chiral extrapolation of this lat¬ 
tice data to the physical pseudoscalar masses using the 
formalism presented in the previous sections, we correct 
for the finite simulation lattice volume using the differ¬ 
ence between infinite-volume integrals and finite-volume 
sums for the loop integral expressions in the chiral expan¬ 
sion. The procedure used here follows Ref. |38|. There 
are no free low-energy constants in the finite-volume cor¬ 
rection expressions; they depend only on A, D and C. 
We fix these constants to their SU(6) values, namely 
D + F = gA = 1.27, F = 2/3£i, C = 2D. We 
also set the pion decay constant in the chiral limit to 
fTT = 0.0871 GeV [32] • The magnitude of the finite vol¬ 
ume corrections is in the range 0.1-0.6% of the relevant 
form factors for each transition and set of pseudoscalar 
masses considered. These shifts, while small, are sig¬ 
nificant compared with the percent-level SU(3)-breaking 
effects of interest here. In all cases the hnite volume cor¬ 
rections act to enhance the size of SU(3)-breaking. To 
account for any model-dependence in our estimate of the 
corrections, we allow F, D and C to each vary by 20% 
about their SU(6) values and add the resulting shift in 
the finite volume corrections in quadrature to the statis¬ 
tical uncertainty on the lattice simulation results. These 
additional uncertainties have a magnitude approximately 
25% as large as the finite-volume corrections themselves. 
To examine the effect of these corrections on our final 
results, we perform the entire fit and analysis described 
below to the uncorrected data set; a comparison is shown 
in the results section. The final extrapolated results with 
and without hnite-volume corrections are consistent at 1- 
sigma for all transitions. 


















ml-ml (GeV^) 



FIG. 3. Vector form factors as a function of the SU(3)- 
breaking parameter — The two (blue and red) trajec¬ 
tories on each plot correspond to the two different constant- 
singlet-mass lines (trajectories {1,2} and 3 in Table|^ respec¬ 
tively) on which the simulations lie. The error bands are as 
described in the text. The vertial dashed line denotes the 
physical psendoscalar mass point. 


After correcting the simulation results to infinite vol¬ 
ume, we fit the free low-energy constants bo and bp to the 
lattice simulation results to perform a chiral extrapola¬ 
tion to the physical pseudoscalar masses. For each choice 
of the constants F, D and C, which are varied within 20% 
about their SU(6) values, and the FRR dipole regulator 
mass A, which is allowed to vary in the range 0.6-1 GeV, 
bo and bp are fit to the numerical simulation results. The 
quality of fit is illustrated in Fig. The central lines of 
each band show the fit using the central values of F, D, 
C and A. The width of the shaded regions is determined 
from the statistical lu variation generated by the fit of 
and bp with the central values of F, D, C and A, added 
in quadrature with the shift in central values obtained 
when fitting bp and bp with the range of choices of F, 
D, C and A. For our final results which are presented in 
the next section, we also add an uncertainty calculated 
by allowing the value of the heavy-baryon mass scale Mq 
to vary between the chiral-limit value and the average 
octet baryon mass at the physical point. 

Finally, we comment that at the physical pseudoscalar 


Bi 

B2 

pS\J{3) 

Jl 

(^1//SU{3) _ 100 Q" 

(GeV") 

E 

N 

-1 

-1.5(ll)Eg(2) : 

-0.065 

H 

E 

1 

-4.8(7) _+5,(4) 

-0.016 

A 

P 

fs 

V 2 

-4.5(7) _+fi3(3) 

-0.031 


A 


-5.4(7)111(0) 

-0.041 


TABLE II. Results for fi{Q^ = —{Mbi — Mb^Y) at the phys¬ 
ical pseudoscalar masses and extrapolated to infinite volume. 
The fourth column shows the SU(3)-breaking corrections as 
a percentage. The first uncertainty is statistical and also in¬ 
cludes the uncertainty in the finite-volume corrections. The 
second uncertainty allows for variation of the low-energy con¬ 
stants D, F and C and the FRR dipole regulator mass A, and 
the third the variation of Mq, as described in the text. 



FIG. 4. Percentage SU(3)-breaking in fi at Q^ax = —{Mbi — 
Mb 2 )^ for each of the strangeness-changing transitions consid¬ 
ered here. The three sets of data points correspond to our full 
results (blue solid squares), as well as the results applying an 
identical analysis but omitting the finite volume corrections 
(orange open circles) or omitting decuplet intermediate states 
(green solid triangles). 


masses we do find nonanalytic corrections to the SU(3)- 
breaking beyond the Ademollo-Gatto terms quadratic in 
the strange-nonstrange quark mass difference. These are 
a consequence of the opening of the meson-baryon decay 
channels as the physical degree of SU(3)-breaking is ap¬ 
proached. For instance, the S —> Al^ amplitude exhibits 
a cusp at the point where the imaginary part associated 
with the decay E —> Ntt appears. The numerical in¬ 
fluence of these terms beyond the standard {nis — mi)^ 
effects tend to be much smaller than the statistical un¬ 
certainties of this calculation. 


V. FIGURES AND RESULTS 

The results for the vector transition form factors, af¬ 
ter infinite-volume and chiral extrapolation, are sum¬ 
marized in Table m Recall that these results are at 
_ Mb^)^ (i-e., corresponding to g = 0 in 
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FIG. 5. Comparison of percentage SU(3)-breaking in /i determined in this work, highlighted by the shaded bands, with that 
of other calculations. The error bands for our results are those given in Table |TT] combined in quadrature. Blue squares, green 
circles and orange triangles denote results of quark model |10llll |. 1/Nc expansion [1^ and chiral perturbation theory |14II16[|40] 
approaches respectively, while the pink diamonds show results from lattice QCD muD]. The red stars show the results of 
this work at = 0 (solid line), where we have corrected from ^ = 0 to = 0 using the dipole form given in Eq. ( |29[ ), and at 
= —{Mb — Mb')^ (dotted line). 


our lattice simulations with fixed zero sink momentum), 
with the physical values of the baryon masses Bi and B 2 , 
instead of at = 0 as is standard. Moving to = 0 
would reduce the magnitude of each form factor, i.e., in¬ 
crease the SU(3)-breaking effect in each case (as will be 
shown explicitly later). As described in the previous sec¬ 
tion, the quoted uncertainties allow for 20% variation of 
the low-energy constants D, F and C from their SU(6) 
values, and for the FRR dipole regulator mass A to vary 
in the range 0.6-1 GeV. Furthermore, we allow Mg, the 
heavy-baryon mass scale used to account for leading rel¬ 
ativistic (or kinematic) corrections, to vary between the 
chiral-limit value and the average octet baryon mass at 
the physical point. We also account for uncertainties in 
the finite-volume corrections as described in the previous 
section. 

Figure 1^ shows the results from Table [Tl] graphically, 
alongside the results obtained using an identical analysis 
but omitting either finite-volume corrections or contribu¬ 
tions from decuplet baryon intermediate states. Clearly, 
all results are stable under these changes. Previous pure- 
effective-field-theory calculations of these quantities (e.g.. 
Ref. [TS]) have typically been very sensitive to decuplet 
baryon effects. We attribute this difference primarily to 
our use of the FRR scheme. 

Following the work in Refs. [551 SIl lll]j we are also 
able to use the chiral extrapolation formalism to deter¬ 
mine the effect of a non-zero light quark mass difference 
(to„ yf rrid) on our results. As we find such charge- 
symmetry violating effects to be one to two orders of 


magnitude smaller than the SU(3)-breaking effects, we 
neglect these differences. Explicitly, we find the differ¬ 
ence in the quantity (/i//f^^^^ — 1) x 100 for —>■ n 

and —>■ p and also 5° —)■ S+ and 5“ —)■ to be in the 
range 0.03-0.04, which is an order of magnitude smaller 
than the statistical uncertainties of our calculation. 

Finally, to estimate the magnitude of the effect caused 
by the non-zero values of used in our analysis, we 
have corrected from = —{Mbi — to = 0 

using the standard dipole parameterisation which is used 
to fit experimental results |43j : 


h{QY 


/i(0) 

(1 + QVM2)2> 


(29) 


where My = 0.97 GeV is chosen, generally universally 
across the baryon octet, for strangeness-changing (and 
0.84 GeV for strangeness-conserving) decays [H]. These 
numbers may be more directly compared with the results 
of previous analyses as shown in F ig.[5l It is clear that the 
naive extrapolation in by Eq. (291 causes a significant 
enhancement of the SU(3)-breaking in our results, partic¬ 
ularly for the S —> V transition where in our calculation 
the value of is the largest. We emphasize that our 
numerical results are presented in Table jll] and obtained 
at non-zero values of Q^; the = 0 results are merely 
shown to facilitate comparison with other work and are 


obtained using Eq. (29) with no attempt to quantify the 
model-dependence of the extrapolation. 

It is clear from Fig. that quark models in general 
predict negative corrections from SU (3)-breaking [THl [U 
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to the vector form factors. This agrees in sign with 
our results although we predict more significant SU(3)- 
breaking effects in all channels, particularly after ex¬ 
trapolation to = 0. Previous analyses using chiral 
perturbation theory [SHH] have typically found SU(3)- 
breaking effects of the opposite sign. We reiterate that 
our analysis differs from these works not only in our use 
of FRR, to which we attribute the stability of our results 
under the inclusion of decuplet degrees of freedom, but in 
additional terms which arise at non-zero values of Q^. We 
interpret these terms, with the free low-energy constants 
fit to the lattice simulation results, as having partially ab¬ 
sorbed the contributions from higher-order effects which 
are omitted from the conventional chiral expansion. 

By providing a framework for the analysis of system¬ 
atic effects in lattice simulations of the vector form fac¬ 
tor, at the simulation values of Q^, this work constitutes 
progress towards a more precise determination of this 
quantity from hyperon semileptonic decays. In turn, such 
a precise determination of /i could lead to an improved 
determination of the CKM matrix element \Vua \ indepen¬ 
dent of extractions from kaon and tau decays. 
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Appendix A: Tables of coefficients 

In this appendix we tabulate the various chiral coeffi¬ 
cients which appear in the expressions of section [l^ 
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